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We study the SU(2) gauge theory on the lattice with different numbers of fermions in the fun- 
damental representation of the gauge group to explore the gauge theory phase diagram. We find 
evidence for an infrared fixed point for ten flavors. The theory with six flavors shows behaviors 
compatible with the existence of a (quasi) stable fixed point, but the large errors in the present 
data do not allow for decisive confirmation of this. 
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1. Introduction 

There exists a class of gauge theories where, under the renormalization group evolution, the 
coupling shows asymptotic freedom at small distances, analogously to QCD, but flows to a fixed 
point at large distances where the theory hence looks conformal. Such theories have applications in 
beyond Standard Model model building. These include unparticles, i.e. an infrared conformal sec- 
tor coupled weakly to the Standard Model and (extended) technicolor scenarios, that explain 
the masses of the Standard Model gauge bosons and fermions via strong coupling gauge theory dy- 
namics [Q, [3| |J. In addition to direct applications to particle phenomenology, the phase diagrams 
of gauge theories, as a function of the number of colours, N, flavours Nf and fermion represen- 
tations, are interesting from the purely theoretical viewpoint of understanding the nonperturbative 
gauge theory dynamics from first principles. While several semianalyitc methods to estimate the 
vacuum phase diagram of a gauge theory exist, the only truly first principle method is constituted 
by lattice simulations. Several initial studies have appeared in literature: for example SU(2) with 
fundamental representation fermions SU(2) with adjoint fermions [g, 0, g, [llj] and SU(3) 
with fermions in the fundamental Jl2| , |i3[ 14] or in the two-index symmetric jl5|], i.e. the sextet, 



representation. 

The studies with Wilson fermions are subject to lattice artifacts proportional to the lattice 
spacing a, and a program to cancel these lattice artifacts has been devised [[If], [l7[ ]. As a motivation 
for this improvement, consider the measurement of the running coupling using the Schrodinger 
functional method: The coupling is measured using a background field and the scale is set by the 
finite size of the lattice. We consider a lattice of volume V = L 4 = (Na) 4 . The spatial links at the 
t = and t = L boundaries are fixed to constant values, while the spatial boundary conditions are 
periodic. The fermion fields are set to vanish at the t = and t = L boundaries and have twisted 
periodic boundary conditions in spatial directions: \y(x + Li) = exp(in/5)\ff(x). At the classical 
level, the boundary conditions generate a constant chromoelectric field and the derivative of the 
action with respect to v\ can be easily calculated: 

dS d - k 

^T = -' (L1) 

on go 

where k is a function of N = L/a and 7} .At the full quantum level the coupling is defined by 



dS\ k 



The perturbative step scaling function defined using the evolution of the renormalized coupling g 
from scale L to scale sL, i.e. 

L(u,s,L/a) = g 2 {go,sL/a)\ g 2 {goX/a)=u 

= u + (-L lfi + -L u N f )u 2 . (1.3) 

The second line gives the formula in perturbation theory to one loop order, and the fermion con- 
tribution is denoted by Ei i. To evaluate these perturbative contributions we use the methods in 



[18, 19], and choose s = 2. The continuum limit of Ei i is given by the fermionic contribution to 
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Figure 1: Contribution of a massless Wilson quark to the step scaling function normalized to its continuum 
value at one loop order in perturbation theory. The top three curves show the result for gauge groups SU(4), 
SU(3) and SU(2) (from top to bottom) for unimproved Wilson fermions, while the lower three curves show 
the result after &{a) improvement has been taken into account. 



the one loop coefficient bo = j8o/(167T 2 ) of the beta function, i.e. 

8= lim Eu/(2MA,iln2), 



(1.4) 



where ^o,i = l/(24;r 2 ). 

The results for the one loop fermion contribution is shown in figure [l] both for unimproved 
Wilson fermions and with ff(a) improvement. One immediately observes that without improve- 
ment, E^i depends strongly on L/a and approaches the continuum limit only for large lattices, 
while with improvement the large lattice artefacts are absent. Clearly this motivates the need to use 
improved actions in the lattice studies of these theories with Wilson fermions. 

In our nonperturbative study we use Symanzik improved Wilson fermions to remedy the dis- 
cretization errors. We measure the running coupling in both models and the mass anomalous di- 



mension in the model with six fundamental fermions. See also [2C] 



2. The model and theoretical tools 



In this section we will introduce the model and some theoretical tools regarding the step scaling 
function and anomalous dimension of the mass operator. We use the basic Wilson Lattice action 

S = S G + S F , (2.1) 

where So is the standard Wilson plaquette action and and Sp is the clover improved Wilson action 



a=l x 



\j/ a (x) (iD + m )y a (x) + ac sw y a (x) - a M V F^ v (x) y a (x) 



(2.2) 



where D is the standard Wilson-Dirac derivative operator including the doubler term. We set the 
improvement coefficient c sw to the perturbative value [21] c sw = 1 +0.1551 (l)go + ^0>o)- ^ e have 
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performed a few short measurements that imply that at strong coupling this is close to the correct 
nonperturbative value for c sw with 6 and 10 fermion flavours. We also include the perturbative 



improvement at the Schrodinger functional boundaries as described in [ |22[ ] 



There is a relation between step scaling funtion ( |1.3| ) and the /3 function: 



-2\n(s) 



a(u,s) j x 



Near the fixed point /3 -function is small and ( |2.3| ) can be approximated with 



21n(2) 



(2.3) 



(2.4) 



We also measure the mass anomalous dimension y = d\nm q / d\n ji of the theory with 6 fermion 
flavours using the pseudoscalar density renormalization constant which is defined as 



Z P (L) 



V3fi 
fp(L/2Y 



where 



1 



12L 6 

fp(xo) = 



12L 6 



d\d\d'yd\.{l'(u)Y,X%'(v)l(y)y 5 X a i;(z)), 

1 r d 3 yd 3 zMxo)r^ a Y(xoK(y)y5h a az)), 



(2.5) 



(2.6) 



(2.7) 



are correlation functions of the pseudoscalar density. Here sources £ and £' are located at the t = 
and t = L boundaries, respectively. For these measurements the boundary matrices at t = and 
t = L are set to unity. The mass step scaling function is then defined as: 



Lp(u,s,L/a) 



sL/a) 



?o,L/a) 



g 2 {go,L/a)=u 



Op{u,s) = lim JLp(u,s,L/a), 



(2.8) 
(2.9) 



and we choose again s = 2. We find the continuum step scaling function G> by measuring Lp at 
L/a = 6 and 10, and doing a quadratic extrapolation. It can be related to the anomalous dimension 
of the mass operator by 



G P {u,s) 



a(u,s) 



do/(2b ) 



exp 



( y(x) d 
^ \P(x) box 



(2.10) 



where bo = j8o/(167T 2 ) in terms of the one-loop coefficient j3o of the beta function and do = 
8/(1 6tt 2 ) is the corresponding one-loop coefficient for the anomalous dimension, y = —dog 2 ■ This 
can be approximated at the fixed point with y*(g 2 ) = — T^j^jjr • 



3. Measurements and results 

The zero mass limit is determined by measuring the fQ. = 1/(8 + 2mo jC ) for all the used values 
of (5 via the PCAC relation using lattice size 16 4 . The measured values of are then used for all 
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lattice sizes. The running coupling g is measured for j8 = {8,6,4,3,2, 1.7, 1.5, 1.3, 1} for Nf = 10 
and p = {8,5,4,3,2.4,2, 1.8, 1.5, 1.44, 1.4, 1.39} foriV/ = 6, using lattice sizes 6 4 , 8 4 , 12 4 and 16 4 . 
The measured values of g are then used to find an interpolating function in j8 of the form 

1 _ f3 " /wv 

gHP,L/a)~2Nto Ci \T) ( } 

with Co = 1. The interpolating function is used to find the step scaling function for L/a = 6,8, and 
the continuum limit is extracted using 

X(u,2,L/a) = o(u,2)+c(L/a) 2 . (3.2) 

Because of the improved action we expect the 6 {a) terms to be subleading. Unfortunately, with 
only two points in the extrapolation, it is not possible to verify the accuracy of the extrapolation 
quantitatively. 

The anomalous dimension of the mass operator is determined similarly: The pseudoscalar 
density renormalization constant Zp is measured for j3 = {2.4,2, 1.5, 1.44, 1.4, 1.39} and for L/a = 
6, 8, 10, 12, 16, 20, and these values are used to find an interpolating function of the form 

Z P (/U/«) = (i) (3.3) 

with co = 1. This is then used to calculate L P (u,2,L/a), which is in turn extrapolated to the 
continuum limit with an order (a/L) 2 ansatz. 

In figure || we show the step scaling functions. In the Nf = 10 theory the evolution of the 
coupling is extremely slow, and our results basically agree with this at g 2 < 2.5: the step scaling 
practically vanishes in this range. In this case we expect the two-loop perturbative step scaling func- 
tion to be fairly accurate, and from figure |2] we see that the errors should be an order of magnitude 
smaller in order to resolve it. At stronger coupling the measured step scaling deviates significantly 
from zero to negative values. Combined with the analytically known weak coupling behaviour, this 
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Figure 3: Anomalous dimension of the mass operator for Nf = 6 

indicates that the /3 -function must have a fixed point somewhere in this range. However, we believe 
that a large fraction of the observed deviation from the perturbative step scaling at strong coupling 
arises from the results at our strongest lattice coupling = 4/gg = 1. This point deviates clearly 
from the rest of the simulation points, possibly indicating stronger cutoff effects. 

In the Nf = 6 theory the evolution of the coupling remains slow, which leads to large errors 
in the step scaling function. In this case we were able to reach mesured couplings up to g 2 = 14. 
However, the results indicate that the possible infrared fixed point is at g 2 > 13, and our statistical 
resolution is not sufficient to confirm or exclude the existance of an IRFP. 

The measured values of mass anomalous exponent y are shown in figure ||. It shows that 
0.1 < 7 < 0.3. Because the value of the anomalous dimension of the mass operator is only scheme 
independent at the fixed point, these results are not that interesting, since we were unable to find 
the IRFP. 

4. Conclusions 

Our simulations verify that the SU(2) gauge theory with 6 flavours of fundamental represen- 
tation fermions is indeed close to the lower edge of the conformal window. Unfortunately, the 
possible fixed point in this theory is at such a strong coupling that we were not able to fully resolve 
the behaviour: the results are compatible either with a fixed point at g 2 > 12 or with a "walking" 
behaviour where the /3 -function almost vanishes. The value of the fixed point coupling is naturally 
scheme dependent; this value is for Schrodinger functional scheme. To resolve this question re- 
quires simulations with an action which can be used at stronger lattice couplings than used in this 
work. 
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